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ABSTRACT 

Context. Massive stars bum hydrogen through the CNO cycle during most of their evolution. When mixing is efficient, or when mass 
transfer in binary systems happens, chemically processed material is observed at the surface of O and B stars. 

Aims. ON stars show stronger lines of nitrogen than morphologically normal counterparts. Whether this corresponds to the presence 
of material processed through the CNO cycle or not is not known. Our goal is to answer this question. 

Methods. We perform a spectroscopic analysis of a sample of ON stars with atmosphere models. We determine the fundamental 
parameters as well as the He, C, N, and O surface abundances. We also measure the projected rotational velocities. We compare the 
properties of the ON stars to those of normal O stars. 

Results. We show that ON stars are usually helium-rich. Their CNO surface abundances are fully consistent with predictions of 
nucleosynthesis. ON stars are more chemically evolved and rotate - on average - faster than normal O stars. Evolutionary models 
including rotation cannot account for the extreme enrichment observed among ON main sequence stars. Some ON stars are members 
of binary systems, but others are single stars as indicated by stable radial velocities. Hence, mass transfer is not a simple explanation 
for the observed chemical properties. 

Conclusions. We conclude that ON stars show extreme chemical enrichment at their surface, consistent with nucleosynthesis through 
the CNO cycle. Its origin is not clear at present. 
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1. Introduction 


Massive stars are born as O and B dwarfs on the main sequence. 
They subsequently evolve into supergiants of various type (blue, 
yellow, red) as their effective temperature decreases and their ra- 
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dius increases. Some stars may evolve back and forth between 
these types of sup ergiants due to as-y et poorly-known physical 
mechanisms (e.g. iGeorgv et al.ll2014ll . Above about 25 Mq (at 
solar metallicity) these stars develop strong winds after the su¬ 
pergiant phase, becom ing nitrogen- or carbon-rich Wolf-Rayet 
stars (WN, WC; e.g., ICrowthen l2007l) . The different types of 
Wolf-Rayet star reflect the different compositions of their sur¬ 
face material: WN stars show nitrogen enrichments and carbon 
depletions, while WC stars are hydrogen-free and have a high 
fraction of carbon in their atmospheres. 

These chemical properties are direct consequences of a dif¬ 
ferent evolutionary states. Massive stars burn hydrogen to helium 
through the partial or complete CNO cycling (depending on the 
temperature). In equilibrium, most of the nuclei in the CNO cy¬ 
cle are in the form of nitrogen; consequently, nucleosynthesis in 
massive stars produces an excess of nitrogen and a depletion of 
carbon and oxygen. If mixing processes are efficient, part of this 
processed material can be brought to the surface, and thus be 
detected spectroscopically. WN stars show the typical chemical 
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patterns of CNO burning; in subsequent evolutionary phases, he¬ 
lium burns to carbon, accounting for the chemical appearance of 
WC stars. 

The CNO cycle proceeds in the cores of massive stars as long 
as hydrogen is available (i.e. during the main sequence). Hence 
we expect massive stars to become increasingly nitrogen rich, 
and carbon/oxygen poor, as they evolve off the zero-age main 
sequence. 

Rotation is a p owerful way of triggering mixing mecha nisms 
in massive stars (iMaeder & Mevneg l2000t iLangerl 12012l) . and 
evolutionary calculations of rotating stars show that the surface 
composition of OB stars can be modified by C NO-processed ma¬ 
terial, even during early evolutionary phases (iBrott et al.ll201 It 
lEkstrom et al.l201^IChieffi & Limongil2013h . The degree of en¬ 
richment depends on several parameters: rotation speed, metal- 
licity, initial mass, magnetic field. 

Observationally, the predictions of stellar-evol ution models 
incorporating rotation have been partly confirmed. iHunter et al.l 
(l2008i 1200^ showed that the majority of the B stars they 
studied in the Galaxy, LMC and SMC exhibit nitrogen en¬ 
richments as predicted by rotating models, although a non- 
negligible fraction (20 to 40%) were fou nd to be more e n riched 
tha n expected for their rotation speed. iPrzvbilla et al.l (1201 Oh 
and iMaeder et al.l (12014h demonstrated that B stars showed sur¬ 
face CNO patterns cons i stent with the expectations of nucleo¬ 
synthesis; iMartins et al.l (12015h reached similar co nclusions for 
a large sam ple of Galactic single O stars (see also iBouret et al.l 

I20lil20nh . 


Surface abundances can also be modified by mass transfer in 
binary systems. If the binary separation is small, the more mas¬ 
sive component fills its Roche lobe first, because of its faster evo¬ 
lution, and may dump processe d material onto the surface of the 
secondary dLanger et al.ll200^ . The primary may subsequently 
explode as supernova, which will either disrupt the system, leav¬ 
ing the chemically-contaminated secondary as a single star, or 
produce a binary system with a high-mass star and a compact 
object. During the mass-transfer process, the primary may also 
transfer angular momentum to the secondary, whic h is thereby 
spun up dWellstein et al.ll^Ol^ iPetrovic et al.ll20d^ . The faster 
rotation can trigger additional mixing of CNO material produced 
in the secondary’s core, contributing further to the modification 
of surface abundances. Even in the absence of mass transfer, the 
rotation of binary components may be affected by tidal interac¬ 
tions, with consequences for mixing. 

Surface abundances are therefore a key to u nderstand i ng the 
evolution o f single and binary massive stars. IWalbornI dl97(]l 
ll97lLfT97^ reported the existence of O and B star s with peculiar 
CNO spectra: the OBN and OBC stars (see also IWalborn et akl 
l2004l) . In the former, lines of nitrogen (especially N iii /i/14630- 
4640) are much stronger than in normal OB stars, while in the 
latter, they are weaker. At the same time, C iii /14650 is weak in 
OBN stars. Most ON stars have a spectral type between 08.5 
and 09.7 - where Cm and Niii lines are easily observed - but 
some are also found at spec tral type 02, based on the m o rphol- 
ogy of N IV and Oiv lines dWalborn et al.l120041) . iLested dl973h 
studied the ON star HD 201345 and concluded that its spec¬ 
troscopic appearance was du e to the presence of CN O pro¬ 
cessed material on its surface. ISchonberner et al.l (Il988h studied 
three ON stars together with two normal O stars, and concluded 
that the ON stars were helium-rich and showed clear signs 
of CNO pr ocessing at their surfac es. Similar conclusions were 
reached by Villatn ^iz et al.l (l2002h for the ON star HD 191423. 
ISmith & HowarthT i 1994 ) determined the helium abundance of 
one OC, one normal O and one ON star. They found an in¬ 


creasing ratio He /H along the OC /O/ON sequence, confirming 
the suggestion of lWalbom (Il976h that the OBN and OBC stars 
represent different degre es of chemical evolution of OB stars. 
iHowarth & SmithI (1200 ih studied the distribution of rotational 
velocities of ON stars and concluded that, on average, they rotate 
faster than normal stars, supporting the hypothesis that rotational 
mixing c ould account for their anomalous surface abundan ces 
(see also ISchonberner et al.lll988h . iBolton & RogersI (Il978h in- 
vestigated the binary frequency of OBN/OBC stars (see also 
iBovaiian et al.l I2 ()(ct . They found that most OBN stars show 
radial-velocity variations (with some binary systems clearly 
identified), while OBC stars seem to be constant. 

While these earlier investigations have provided evi¬ 
dence of CNO processing, to date qu antitative spectro¬ 
scopic analyses have focussed on helium ([Schonb erner et all 
119881: ISmith & Howarthl[l9^ IVillamariz et ^ 20021) . emplov- 
ing relatively simple model atmospheres; determinations of 


CNO abundances have relied on e quivalent-width measure¬ 
ments and curve-of-g rowth analyses (ISchonberner et all 119881: 
IVillamariz et al.l 12002 ). Consequently, the status of OBN/OBC 
stars remains unclear, and the origins of the chemical-abundance 
peculiarities observed at their surfaces are not fully understood - 
e.g., what are the relative roles of rotation and binarity? In order 
to make progress in the understanding of OBN stars, we have 
conducted the first quantitative analysis of the helium and CNO 
abundances of a significant sample of ON stars, using modern 
non-LTE/line-blanketed atmosphere models. Sectionl^of the pa¬ 
per presents the sample and the observational material; analysis 
methods are described in Section[3 and the results in Section|4] 
These results are discussed in Section |5j with our conclusions 
given in Section|6] 


2. Observations and sample 

The most complete catalogu e of Galactic l ate-typ e ON stars 
known to date is provided bv lWalborn et alJ (1201 ll) . It contains 
13 stars: four dwarfs/sub-giants, seven giants, and two super- 
giantsQ In the present study, we consider twelve stars from this 
sample. We exclude the ON9.7 Il-III(n) star HD 89137 because 
it is a double-lined spectroscopic binary and the available data 
are not sufficient to separate unambiguously the two compo¬ 
nents. 

The spectroscopic data were collected from a number of 
archives and unpublished material; sources and dates of obser¬ 
vation are listed in Table [T] We give a brief description of each 
dataset below: 

- AAT data were acquired at the 3.9-m Anglo-Australian Tele¬ 
scope with the University College London Echelle Spectro¬ 
graph (UCLES) and a 31.6 lines mm ' grating. Observations 
were made at several spectrogrograph settings in order to 
obtain near-continuous wavelength coverage ~3650-7000A 
at a resolving power R ^ 23 000, and signal-to-noise levels 
>100 throughout. 

- The lACOB spectroscopic database (ISimon-Dlaz et al.l 
1201 ll) includes multi-epoch observations for seven of the 
stars listed in Table [1] The observation s, obtained with the 
Elbre-fed Echelle Spectrograph (LIES, iTelting et~^l2014l) 

* IWalborn et all ll201 ll) explicitly excl ude two f urther stars. The first, 
BD-l-36 4063, is an 09.7 lab supergiant iMathvsl ( I19891) . known to be a 
spectroscopic binary dWilliams et alJ2 009l). The second, HD 105056, is 
classified as ON9.7 lae (e.g.. lMai^Ap ellaniz^tjilj2004l). but it may be 
a lower-mass, post-AGB object dWalborn et alJ201 llf 
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Table 1. Observational information. 


ST Instrument date of observation 


Star 


HD 12323 
HD 13268 
HD 14633 
HD 48279 
HD 91651 
HD 102415 
HD 117490 
HD 123008 
HD 150574 
HD 191423 
HD 191781 
HD 201345 


ON9.2 V 
ON9.5 Illn 
ON8.5 V 
ON8.5 V 
ON9.5 Illn 
ON9 Ilinn 
ON9.5 Ilinn 
ON9.2 lab 
ON9 Ill(n) 
ON9 Il-IIInn 
ON9.7 lab 
ON9.2 IV 


TIES 

FIES 

ESPaDOnS 

EIES 

EEROS 

EEROS 

EEROS 

HOLES 

EEROS 

ELODIE 

HOLES 

ESPaDOnS 


08 sep 2011 
average 
09 oct 2009 
average 
05 may 2009 
10 jun 2008 
14 may 2008 
22 jun 1992 
10 jun 2008 
29 aug 2004 
14 aug 1995 
26jul 2010 


Notes. Spectral types are from ISota et all d201 iLlTOlih . 


attached to the 2.5 m Nordic Optical Telescope (NOT), cover 
the spectral range 3700-7300A with R=46 000. Data reduc¬ 
tion was performed with FIEStooQ and the normalization of 
the spectra using own procedures developed in IDE. 

- Spectra of four targets were extracted from the ESO/FEROS 
archive. These data were ob t ained in the context of the OWN 
project (iBarba et al.l 1201 Ol l2014ll . They have a resolving 
power of 48 000 and were reduced by the fully automated 
pipeline distributed by ESO. 

- HD 191423 was observed with the Elodie Spectrograph, 
mounted on the 1.93-m telescope at Observatoire de Haute- 
Provence (OHP). The resolving power is 42000 and the 
wavelength coverage 3895-6815 A. Exposure times were 
chosen to ensure a signal-to-noise ratio of at least 100 at 
5200 A. Reduction was p erformed using the st andard reduc¬ 
tion pipeline described in iBaranne et al.l (Il996h . 

- For HD 14633 and HD 201345, spectra were acquired with 
the ESPaDOnS spectropolarimeter mounted on the Canada- 
France-Hawaii Telescope. The resolving power is 65 000 
and the wavelength coverage is 3700-10500 A. Data reduc¬ 
tion was performed with the automated procedure Libre Es¬ 
prit, a full description of the data is given by Wade et al. 
(2015, submitted). 

Table[T]lists the observational data used for the spectroscopic 
analysis. Additional data are described in Table lA.il and will be 
discussed in Sect. 15.31 

3. Modelling and spectroscopic analysis 

We have used the atmosphere code CMFGEN to analyze the sur¬ 
face properties of the ON stars. A full description of CMFGEN 
can be found in iHillier & Mili^ (Ho^ . In a nutshell, CMFGEN 
solves the radiative transfer and statistical equilibrium equations 
in the comoving frame, leading to non-LTE models. The temper¬ 
ature structure is set from the radiative-equilibrium constraint. 
Spherical geometry is adopted to take into account extension due 
to the strong winds of O stars. The density structure is computed 
from mass conservation and the velocity structure is constructed 
from a pseudo-photospheric structure connected to a jS-velocity 
law. The photospheric structure is obtained from a few iterations 

^ http://www.not.iac.es/instruments/fies/fiestool/ 


of the hydrodynamical solution in which the radiative force com¬ 
puted from the level populations and atomic data is included. The 
■final synthetic spectrum is obtained from a formal solution of the 
radiative transfer equation. 

The synthetic spectra are subsequently compared to obser¬ 
vations to determine the stellar parameters. In practice, we pro¬ 
ceeded as follows: 

- Rotation and m acro turbulen ce: we used the Fourier- 
transform method (lGravlll97^ to determine V sin/, the pro¬ 
jected equatorial rotational velocity. We relied on O m 45592 
when possible; otherwise, we used Hei 44920. We subse¬ 
quently estimated the expected and log.? from the tar¬ 
get’s spectral type, using the calibration of iMartins et al.l 
(l2005h . and selected a corresponding synthetic spectrum 
form our database of models. We convolved this spectrum 
with a rotational profile (adopting V sin/ from the previ¬ 
ous step), and performed an additional convolution by a 
radial-tangential profile (parameteri zed by a velocity Vmac) to 
take m acroturbulence into account dSimon-Dlaz & Herrerol 
120141) . We varied Vmac until a good match was obtained to 
the observed spectrum (especially around 4100 A, 4700 A 
and 4900 A). Hncertainties on V sin/ and Vmac are of the or¬ 
der ~10 and 20 km s ' respectively. 

- Effective temperature: we relied on the traditional ionization- 
balance method to constrain Teff. We used the helium lines 
He 1 44026, He ii 44200, He 1 44388, He 1 44471, He ii 44542, 
He I 44712, Hei 44920, Hei 45016, and Hen 45412. 
We found that when Hen 44542 was perfectly matched. 
Hen 44200 was usually slightly too strong and Hen 45412 
slightly too weak. This is partly attributed to the echelle na¬ 
ture of most of our spectra and thus to uncertain normal¬ 
ization. As a consequence, effective temperatures are deter¬ 
mined within about +1500 K. 

- Surface gravity: the wings of Balmer lines were used to de¬ 
termine log g, with a typical uncertainty of ~0.15 dex. 

Since our principal focus is on surface parameters, we did not 
try to reproduce the details of emission lines likely to be formed 
in the wind (particularly Hen 44686, Ha); instead, we simply 
adopted combinations of terminal velocity and mass-loss rate 
that give a reasonable fit to these features in each star, with the 
wind-acceleration parameter f fixed at 1.0. Luminosities were 
adopted from the calibration given bv IMartins et al.l (l2005l) . In¬ 
strumental broadening was negligible given the high resolution 
of the observed spectra and the generally large rotational broad¬ 
ening, and our final synthetic spectra were simply convolved by 
rotational and radial-tangential profiles (see above). 

For the determ i nation of surface abundances, we proceeded 
as in iMartins et al.l (l2015h . For a given Teff and log g we ran mod¬ 
els with different He, C, N, and O abundances. We identified the 
cleanest lines of each element in the observed spectrum, and per¬ 
formed analysis combining all selected features. The com¬ 
puted values were renormalized to a minimum value of 1.0, 
with the best-fit abundance taken to correspond to this minimum, 
and 1-cr uncertainties defined by;^f^=2.0. 

For helium, we utilized from seven to ten lines selected from 
Hei 44026, Hen 44200, Hei 44388, Hei 44471, Hen 44542, 
He 1 44713, He 1 44920, He 1 45016, He ii 45412, and He 1 46680. 
Figure[T]shows an example of He/H determination, illustrated by 
HD 48279. An abundance greater than 0.1 (by number) is clearly 
favoured, simply by inspection. The quantitative analysis restab- 
lishes that He/H=0.16+0.05 best reproduces the set of helium 
lines. 
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Fig. 1. Determination of the ratio He/H for HD 48279. Observed He i 
and Hen lines are shown by the black solid line. The blue dotted 
(red dashed; green dot-dashed) lines are models with He/H=0.10 (0.15; 
0.30). The lower right panel shows the^^ results from the quantitative 
analysis, from which one determines He/H=0.16±0.05. 


Our primary carbon-abundance diagnostic was C iii AAQIQ. 
For a few stars. Cm /14163, Cm /14188, and Cii 44267 could 
also be used. We did not use C iii 44650 because of its sensi¬ 
tivity to_wmdvto_memllic^ and to uncertainties in the atomic 
data (iMartins & Hilliei|l2012h . In most cases, we could obtain 
only an upper limit on C/H because C iii 44070 is the only use¬ 
ful carbon line that could be detected, and is rather weak. Dif¬ 
ficulties in normalizing the spectra resulting from the presence 
of nearby Hb further undermine any attempt to put a better con¬ 
straint on C/H. Figurel^shows an example of C/H determination, 
for HD 117490. In this case a conservative upper limit on C/H 
of 3 X 10“^ is adopted. 


Fig. 2. Observed C iii 4070 line of HD 117490 (black solid line) together 
with models with different C/H values (colored interrupted lines). 

We adopted a photospheric microturbulent velocity of 
10 km s ' in our synthetic spectral A larger value 
would ten d to reduce the der i ved abundances, especially 
of helium (iMcErlean et alJ Il998l : IVillamariz & Herr^ l200(i 
lHowarth& Smithll2001h . Using HD 191423 as a test case, we 
find that increasing the microturbulence to 20 km s ' leads to a 
reduction of 0.05 in the inferred He/H ratio. 

As a check on our results, we performed a separate analysis 
of ou r sample using the iacob-gba t package (ISimon-Dfaz et al.l 
1201 ll : l^bln-Saniulian et alJl2014l) . which uses a -fitting al¬ 
gorithm coupled to a large, pre-computed grid of F ASTWIND 
models (ISantolava-Rev et^ ll997tlPuls et al.ll2005h . This com¬ 
pletely independent approach allows the determination of Teff, 
\ogg, V sin/, Vtnac, and He/H. The results turned out to be in 
excellent agreement with the CMFGEN analyses, within the re¬ 
spective error bars, lending strong support to our quantitative re¬ 
sults. 


The principal nitrogen-abundance indicator was the N iii 
4510^535 A complex of lines. Occasionally, other indicators 
were added: Nii 43995, Nii 44004, Niii 44044, Nii 44447, 
Nil 44607, Nil 45001, Nii 45005, Nii 45011, Nii 45676 and 
NII45680. Finally, O iii 45592 was adopted as the main oxygen- 
abundance indicator. In some stars, we could also make use of 
On 43913, On 43963, On 44277-78, On 44318, On 44368, 
On 44416-18, On 44603, and On 44611. Compared to normal 
O stars, fewer useful lines are generally available to determine 
surface abundances beca use the ON stars typica lly have rela¬ 
tively large V sin/ values (iHowarth & Smithll2001h : hence weak 
metallic lines are too broad to be readily detected. Lines from 
two ionization states are usually taken into account in the abun¬ 
dance determination, so that the error bars include uncertainties 
related to flawed N and O ionization balances. At the same time, 
by performing the abundance determination on an ensemble of 
lines we minimize uncertainties related to the forma tion/physics 
of individual lines (see Fig. 1 in iMartins et al.ll^l5h . 


4. Results 

The results of the spectroscopic analysis are summarized in Ta¬ 
bled and the best-fit models compared to the observed spectra in 
Figs. IB.lllB.12l We emphasize that the numerical results listed 
in Table |2] should be viewed as ‘surface average’ values; many 
of the sample stars have large rotational velocities, which can 
affect their shap es, and create surface gradients in the phy sical 
parameters (e.g. lHowarth & Smithll200ll : lPalate et alJl201^ . For 
HD 191423, the fastest rotator of the sample, rotational broaden¬ 
ing is so large that normalization of the echelle spectra was ex¬ 
tremely difficult, and carbon and oxygen abundances could not 
be determined. 

Table El shows tha t most ON stars are helium-rich. 
ISmith & HowarthI d 19941) performed a helium-abundance deter¬ 
mination for HD 123008 and obtained He/H= 0.20 + 0.05, in 

^ Except for HD 123008 for which we found that 20 km s“‘ gave a 
better fit of the helium lines 
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Table 2. Parameters of the sample stars. 


Star 

Spectral 

Teff 

logg 

loggc 

V sin; 

^mac 

He/H 

C/H 

N/H 


Type 

[kK] 



[kms-'l 

[kms-'l 


[10-4] 

[10-4] 

HD 12323 

ON9.2 V 

33.5 

4.00 

4.01 

130 

- 

006 

<0.3 

3.5!f' 

5 0+“-^ 
-2Q 
5.2//° 

HD 13268 

ON9.5 Illn 

32 

3.50 

3.63 

310 

- 

0.20+0.10 

<0.5 

HD 14633 

ON8.5 V 

34 

3.80 

3.81 

100 

120 


<0.3 

HD 48279 

ON8.5 V 

34.5 

3.80 

3.82 

137 

50 

0.16+0.05 

U.UJ_0.20 

^■^-2.0 

HD 91651 

ON9.5 Illn 

31 

3.50 

3.62 

310 

- 

0.14/008 

<0.5 

5 4+1.6 

HD 102415 

ON9 Ilinn 

31 

3.50 

3.70 

376 

- 

0.21+0.10 

<0.6 

7.6+11 

HD 117490 

ON9.5 Ilinn 

30.5 

3.50 

3.66 

375 

- 

0 05 

<0.3 

HD 123008 

ON9.2 lab 

30 

3.10 

3.10 

37 

100 

0 21+° '' 
'+■^ + -0.07 

0 22+° “' 

13.5/1^ 

HD 150574 

ON9 Ill(n) 

31 

3.40 

3.49 

240 

- 

0.23+0.06 

<0.5 

>10.0 

HD 191423 

ON9 Il-IIInn 

31.5 

3.50 

3.72 

445 

- 

0.25+0.08 

- 

>5.0 

HD 191781 

ON9.7 lab 

28 

3.10 

3.12 

107 

70 

0-251? 

<0.5 

7 3+''-° 
-5Q 
4 0+^-^ 

-1.3 

HD 201345 

ON9.2 IV 

34 

4.00 

4.01 

95 

60 

0.1 

<0.4 


O/H 
[ 10^1 
2.5+^-^ 

3.1; 

1.7 
4 

-?-3 
9 0 + 1.9 

<3.0 

2 

- 2-0 

4.8£ 

6-0!z3 


U 


o 1 +0.8 


Notes. Uncertainties on Tcff, logg, V sin; and Vmac are ~1.5kK, 0.15 dex, 10 and 20 km s ' respectively, logg^. is the surface gravity corrected for 
centrifugal acceleration. Abundances are number ratios. 



log (U„) 


Fig. 3. logg - logfTeff) diagram for the sample stars. Triangles (pen¬ 
tagons, squares, circles) are for luminosity class V (IV, III, I) stars. Typ¬ 
ical uncertainties are sh own in the upper-left c orner. Evolutionary tracks 
including rotation, from lEkstrdm et akl (l2012h . are overplotted, labelled 
by ZAMS masses. 

good agreement with our estimate. iHowarth & SmithI (1200 ih de¬ 
termined He/H in HD 191423 using simple non-LTE models 
without line-blanketing but a better treatment of 2D effects than 
our approach. They obtained He/H= 0.23 + 0.04 in very good 
agreement with our value. The relatively large helium enrich¬ 
ment of most ON stars indicates a peculiar chemical history 
compared to normal OB stars, whic h do not show system atically 
values of He/H larger than 0.1 (e.g. lMokiem et al]l2007h . 

Figure [3 shows the surface gravity as a function of effec¬ 
tive temperature in our sa mple. ’Geneva’ evolutionary tracks, 
from lEkstrom et aP (l2012l) . are overplotted. We see that the ON 
stars follow a relatively clear sequence; dwarfs have higher logg 



He/H 


Fig. 4. N/H as a function of He/H for the ON stars. The abundance 
ratios are by number. Symbols have the sa me meaning as in Fig . |3 
The evolutionary tracks including rotation of lEkstrbm et all ( 120121) are 
overplotted. The bold part corresponds to the main sequence (central 
hydrogen mass fraction > 0). 

values than giants, while supergiants have the l owest logg, a 
trend similar to that found for normal O stars by iMartins et al.l 
(l2015h . Stars classified as ON have a rather narrow range of ini¬ 
tial masses; 20-25 Mq for the dwarfs/giants, and a little over 
25 Mq for the supergiants, according to the Geneva tracks. 

Figure |4] illustrates the N/H ratio (by number) as a function 
of He/H. The evolutionary models predict a correlation between 
N/H and He/H that simply results from nucleosynthesis through 
the CNO cycle; the higher the helium content, the higher the ni¬ 
trogen abundance. Figure |4] shows a possible trend of this kind 
among the ON stars, but the rather large error bars on the abun¬ 
dance determinations don’t allow for a clear-cut conclusion ob- 


Article number, page 5 of[l5] 
























































A&A proofs: manuscript no. on 



Fig. 5. The N/C abundance ratio as a function of surface gravity for 
the ON stars (large symbols), compared to results for morphologically 
normal stars with sp ectral types 08.5-9.7 (small symbols), taken from 
iMartins et al.l l l201^ . 

servationally. Most stars appear more enriched than expected 
from their position in the logg - log(7'eff) diagram. 

Figure |5] shows the ratio of nitrogen to carbon surface abun¬ 
dance as a function of surface gravity for the ON star^ and 
for a reference sample of morphologically normal O stars taken 
from IMartins et al] (l2015h . For the reference sample, only stars 
with spectral types between 08.5 and 09.7 are considered, cor¬ 
responding to the spectral-type range of the ON sample. The 
dwarf, giant, and supergiant stars are well separated in this plane, 
among both ON and reference groups, largely reflecting the 
surface-gravity differences between luminosity classes, as is also 
seen in Fig. |3] The most striking feature of Fig. |5]is the higher 
log(N/C) ratios found for ON stars, for a given luminosity class. 
The difference amounts to 0.5-1.0 dex or more, since most of the 
log(N/C) values for ON stars are actually lower limits (due to 
the upper limits on C/H; Table|2]i. Another possible trend is that 
of higher log(N/C) among ON stars as one moves from dwarfs 
to giants and supergiants. This tr end was clearly estab lished for 
morphologically normal stars by IMartins et al.l (1201 5h . and ap¬ 
pears to hold also for ON stars. This suggests that the physics 
of chemical mixing follows the same patterns in normal and ON 
stars, pointing to a common origin for surface chemical enrich¬ 
ment. However, we caution that the trend for ON stars is less 
clear due to both the small number of objects, and the fact that 
for most stars we have only lower limits on log(N/C). 

Figure |6] shows the log(N/C) versus log(N/0) diagram. The 
upper-left panel reveals that the ON stars are syste matically more 
enrich ed than the comparison O stars taken from IMartins et al] 
(l2015l) . They extend the relation between log(N/C) and log(N/0) 
observed for normal stars to higher abundance ratios, and are 
also located between the limits corresponding to partial CN 
burning and complete CNO burning. Consequently, the abun- 

As discussed previously, only upper limits to C/H are available in 
most cases. 


dance patterns of ON stars are consistent with the predictions of 
nucleosynthesis through the CNO cycle. 

Other panels of Fig. |6] compare results for ON and normal 
O stars at different luminosity classes. In each case, the ON 
stars are clearly separated from the comparison stars, indicating 
that their surface abundances result from a much stronger mix- 
ing than that experienced by normal stars. IMartins et al.l (12015h 
showed that, on average, supergiants are more enriched than gi¬ 
ants and dwarfs. We saw in Fig.|5]that this trend may exist among 
ON stars, too. Fig. |6] tends to confirm that there is a sequence of 
higher enrichment when moving from dwarfs to giants and su¬ 
pergiants, although there is one outlier (HD 14633, a dwarf that 
appears to be as chemically mixed as the giants). 

Figures |5] and |6] show gaps between the distributions of the 
ON stars and those of the comparison stars, and the reader may 
wonder whether these are real features. In this context, then 
supposing that chemica l mixing may be due to rotation (and 
iHowarth & Smithll2001l have shown that ON stars rotate on aver¬ 
age faster than normal O stars), it is informative to consider the 
distributions of projected rotational velocities in the two sam¬ 
ples. This comparison is performed in Fig. |2l clearly, on aver¬ 
age the stars in the comparison sample are rotating much more 
slowly than those of the ON stars. The median V sin/ value for 
the ON and comparison-O samples are 208.0 and 52.5 km s ', 
and a KS test indicates a probability of less than 1% that the two 
populations are drawn from the same parent distribution. 

Fig. |7] further includes the V sin/ distributio n of the lACOB 
sample studied bv ISimon-Diaz & Herrerol(l2014l) . They provided 
projected rotational velocities for 199 Galactic O stars, which 
can be viewed as a reference distribution for O stars in general; 
the ON stars have highe r V sin/ values, and our comparison O 
stars lower ones, than the lSimon-Dfaz & Herrerol reference sam¬ 
ple. Potentially, therefore, the very different ranges of V sin/ cov¬ 
ered by our two samples (ON and comparison) can explain the 
gap seen in Fig. |6l if rotation leads to a continuous increase of 
mixing as V sin/ increases, and if our two samples are well sep¬ 
arated in terms of projected rotational velocity, then a separation 
into two groups in the log(N/C)-log(N/0) and log(N/C)-log g 
diagrams is expected. 

5. Discussion 

5.1. Spectral classification and ON stars 

The spectroscopic classification qualifier ‘N’ was developed in 
the context of late-O/early-B stars, as is illustrated by our sam- 
ple, which is l i mited to spectral types between 08.5 and 09.7. 
IWalborn et al.l (l2004l) discovered a related morphological di¬ 
chotomy in terms of CNO features at the earliest spectral type, 
02, where some stars show relatively strong O iv lines and weak 
N IV lines. The opposite behaviour is found in a second group, 
for which the spectral classification ON2 was therefore created, 
by analogy to the behaviour in late-O stars. 

It is important to note that different lines are used when clas¬ 
sifying ON2 and late-type ON stars; lines from more highly ion¬ 
ized elements are observed in the former group (O iv, N iv) than 
in the latter (C iii, N iii). This, rather than some fundamental dif¬ 
ference in surface chemistry, largely accounts for why the ‘N’ 
qualifier is not assigned to O stars at intermediate spectral types 
- when moving from 02 to later spectral types, the diagnostic 
Oiv and Niv lines disappear because of ionization effects.On 
the other side of the classification scheme, the C iii and N iii 
lines around 4630^650A, which define the late-type ON cate¬ 
gory, progressively go into emission at spectral types earlier than 


Article number, page 6 of[l5] 


















































Martins et al.: Abundances of ON stars 





Fig. 6. Distributions of log (N/C) versu s log(N/0) for the O N stars (large symbols) and for a comparison set of normal O stars with spectral types 
08.5-09.7 (small symbols) taken from iMartins et al] ( l2015h . The upper-left panel shows the full sample. The upper-right (lower-left, lower-right) 
shows the dwarf-subgiant (giant, supergiant) sample. Solid lines in the upper-left panel show the relations expected for partial CN or complete 
CNO burning in equilibrium. 


08. In principle, the strength of these emission features could 
be used in a similar manner to their absorption counterparts in 
late-type stars, but the strength o f the emission in these lines 
is complicated by win d effects (iRivero Gonzalez et al. 201 lb 
iMartins & Hilliedl2012l) . Interestingly. IWalborn et al.l (l2010h re- 
cently dehned the Ofc class, corresponding to stars showing 
N III /14630-4640 in emission and C iii /14650 of comparable 
strength. The difference between Ofc and Of stars may be the 
equivalent of the ON phenomenon at 02 and late-O spectral 
types (Of stars having comparatively weaker C lines than Ofc 
stars). 


5.2. Rotation and the ON phenomenon 

Fast rotation could explain the strong chemical enrichments 
that we have found in O N stars, which rotate on average 
faster than normal O stars (iHowarth & Smith! 1200 iL and sec¬ 
tion 00 Theoretical predic tions indicate that faster rotation 
leads to stronger mixing (iMaeder & Mevnet 200r)t Langei 
201^ ; iBrotteyik 2()1 U ^strom et al.ll2012 ’ 1chieffi & Limongi 


20131 : iGeorgy et alJ I2013al) . and thus to the appearance of 
more-strongly proce ssed material at the surface of OB stars. 
IWalborn et al.l ( 2004h performed surface-abundance determina- 


^ We recall that such fast rotation could result from formation pro¬ 
cesses, or may arise through tidal interactions in binary sy stems (which 
would not necessarily have experienced mass transfer; iLanger et al.l 
l2Q08h . 
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Vsini [km s"‘] 


Fig. 7. Cumulative distribution functions of projected rotational veloc¬ 
ities for the ON sample (red circles), o ur comparison sample (black 
squares), and the lACOB 0-star sample dSimon-Diaz & Herreroll20T3 : 
orange triangles). 


tions for 02 and ON2 stars and concluded that the latter were 
more chemically processed than the former. From their positions 
in the Hertzsprung-Russell diagram, they argued that ON2 stars 
may be the product of homogeneous evolution, which is usu¬ 
ally understo od as a consequence of fast rotation (iMaedeill 19871: 
lLangeriri992h . 

Figure [8] confronts our abundance determinations with sev¬ 
eral stellar-evolution models, comp u ted fo r 25 Mp stars at so¬ 
lar me tallicity, from lEkstrom et al.l (l2012l) . IChieffi & Limongil 
(1201 3h . and a dedic ated model comp uted with the STAREVOL 
code (ISiessI 120061: iDecressin et al.l l2009h for which we im¬ 
plemented the same physical ingredients as in lEkstrom et al.l 
BoH : core overshooting, mass-loss prescriptions, angular- 
momentum-transport equation, and turbulent-transport prescrip¬ 
tions. All the models start on the zero-age main sequence with 
an equatorial surface velocity of 300 km s * (equivalent to 
V 0.414, where 14 the critical surface velocity). 

All the inferred abundances lie close to the theoretical 
tracks, which confirms that they broadly follow the predictions 
of nucleosynthesis and chemical mixing However, the abun¬ 
dances predicted at the end of the main sequence (i.e., when the 
hydrogen mass fraction in the core reaches zero) vary signifi¬ 
cantly between codes, due to the sensitivity of the mixing effi¬ 
ciency on both the a dopted prescrip ti ons for the shea r turbu lence 
(vert ical shear frorn lMaedend 19971) in Ekstrom et al. (|^ 012h . and 
from iTalon & Zahnl(ll997 ) in Chieffi & LimongJ i 2013h ). and on 


the numerical treatment itself. Nevertheless, we see that the ON 
stars are all located beyond the end of the main sequence in the 
log(N/C)-log(N/0) plane. While this may be possible for super¬ 
giants, the giants and dwarfs are certainly expected still to be 
core hydrogen burning (see Eig. |3ll. The surface abundances of 
ON stars are thus more enriched than predicted by any of the 


® The STAREVOL calculations were stopped just after the end of the 
main sequence. 


available models, assuming a standard rotation rate on the main 
sequence. _ 

Scrutiny of 15-Mo models from iGeorgy et al.l (1201 3bl) sug¬ 
gests that surface enrichment of CNO-cycle products can be in¬ 
creased significantly (by up to a factor of 3 compared to normally 
rotating stars) through adopting initial rotation rates close to the 
break-up velocity. It is therefore tempting to speculate that the 
nitrogen abundances of ON stars can be attributed to enhanced 
efficiency of rotational mixing in extremely fast-rotating stars. 
To test this hypothesis, new 25-Mo models were computed us¬ 
ing the Geneva code, with ratios of initial to critical equatorial 
velocities, V/Vc, of 0.6 and 0.8 (see Eig. [3]l. Results are pre¬ 
sented in Eig. 0 The left-hand panel shows that increasing the 
initial rotational velocity leads to a more ‘upward’ evolution 
(i.e., the effective temperature decreases less at higher velocity). 
All non-supergiant ON stars are located between the zero-age 
main sequence and the terminal main sequence, whatever the ini¬ 
tial rotation speed. The middle panel of Eig. |9] shows the corre¬ 
sponding surface abundances; as expected, faster rotation leads 
to stronger enrichment. At the end of the main sequence, the 
V/Vc - 0.8 model barely reaches the region occupied by the ON 
stars. However, this happens when logg < 3.5; for the range of 
surface gravities of non-supergiant ON stars (3.5^.0) the sur¬ 
face enrichment is weaker. Consequently, evolutionary models 
with high rotation velocities do not appear to explain the extreme 
enrichment of ON stars, at least with the current formalism used 
to treat rotation. 

The right-hand panel of Eig. |3 shows the effect of adding 
magnetism to the V/Vc - 0.6 model. This models includes the 
effects of the internal magnetic fields on the mixing of chemicals 
and ang ular moment u m, according to the Tayle r-Spruit dynamo 
theory (ISpruitl 120021: iMaeder & Mevnej l2004l) . In this frame¬ 
work, internal magnetic fields produces a strong coupling be¬ 
tween the core and the surface, and considera bly increases the in¬ 
ternal mixing (e.g. lMaeder & Mevnelil20d5h . Significant enrich¬ 
ment, consistent with that of ON stars, can be produced while 
logg >3.5. However, inspection of the left-hand panel of Eig. 0 
reveals that this model produces a blueward evolution and never 
reaches the position of the ON stars. 

Models aside, if rotational mixing is the principal cause of 
the ON phenomenon, one may wonder why not all fast-rotating 
O stars are classified as ON. Eor example, f Op h has a projected 
rotati onal velocity of about 400 km s ' (e.g., iMarcolino et al. 
2009), but its spectral type is simply 09.2 IVnn dSota et al. 
20141) . where the ‘nn’ qualifier reflects the strong rotational 
bro adening; that is, i/ Oph i s not classified as an ON star (see 
also lHowarth & Snfldill2001l) . IVillamariz & Herrerol (l2005l) per¬ 
formed a quantitative analysis of its stellar parameters and sur¬ 
face abundances (using eq uivalent widths measure ments and the 
curve-of-growth method). IMarcolino et al.l (l2009l) . using more 
recent atmosphere models than Villamariz & Herrero, refined the 
Tcff and log g deter mination. C Oph, toge ther with the fast rota¬ 
tors of the sample of iMartins et al.l(l2015l) (i.e. stars with V sin/ > 
250 km s ') are shown in Eig. [10] The left-hand panel of Eig.lTOl 
demonstrates that the ON stars are more chemically processed 
than the non-ON fast rotators. The only non-ON fast rotator with 
a strong enrichment is HD 192281 (04.5 Vn(f)), the most mas¬ 
sive dwarf of the sample; its enrichment is proba bly due to its 
higher mass (see discussion in iMartins et al.ll^l5l) . 

Since chemical enrichment depends not only on rotation but 
also on metallicity, mass and age, one may wonder if the non- 
ON fast rotators are less evolved or less massive than the ON 
stars (supposing their metallicities to be similar since all stars are 
located in the Galaxy). The right-hand panel of Eig. ITOl shows 
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log (N/0) 



M = 25 Mq 
v=300 km/s 


Main sequence 


Post-main sequence 


Fig. 8. As Fig. [^ with evolutionary tra cks for a ZAMS mass of 25 Mp, inc luding rotation (with an initial equatorial velocity of about 300 km s ‘). 
Models are fro mEkstrbm et al.l l l2012ll (top leftl. lChiefh & Limongil ( l2013h (top right), and the STAREVOL code (bottom left). Thick line sections 
correspond to core hydrogen burning. 


that ON and non-ON stars are distributed in the same area of 
the log.g-reff diagram. Thus, they have similar mass and age 
ranges. In conclusion, if ON stars appear to rotate on average 
faster than normal stars, not all fast rotators necessarily display 
the ON phenomenon. Another mechanism may be required, in 
addition to rotation, to produce ON stars. 

5.3. Binarity and the ON phenomenon 

iBolton & RogersI (Il978h reported that most of the ON stars 
known at the time of their study were radial velocity (RV) vari¬ 


ables, but OC stars were all RV constant. IWalborn et al.l (1201 Ih 
discussed radial-velocity variability and the presence of com¬ 
panions in the most recent ON sample. Of the two O super¬ 
giants in our study, HD 123008 is reported to be constant in 
radial velocity, while HD 191781 may be variable. The vari¬ 
ability status of the giants is quite uncertain, with some claims 
of RV variations but based on limited numbers of measure¬ 
ments. The nature of the tentative variability observed in the 
giant sample is also unclear: binarity or wind/photospheric vari¬ 
ability? Of the five dwarfs/subgia nts, HD 48279 is constant in 
radial velocity (iMahv et alJl2009h . HD 12323 is a single-lined 
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Fig. 9. Effect of rotation on a 25 Mq evolutionary model. Left: log g - log(refj) diagram. Middle: log(N/C) - log(N/0) diagram for models at different 
velocities. Right: log(N/C) - log(N/0) diagram for a model including magnetic field and V/Vc=0.6. The bold part of the tracks corresponds to the 
main sequence. Dots in the middle panel correspond to a surface logg = 3.5, i.e. the mi nimum value for dwar fs and giants in the left-hand panel; 
in the right-hand panel, logg is always higher than 3.5. Models for V/Vc = 0.4 are from lEkstrdm et aP i2012h . Models for higher velocities have 
been computed for the present study with the Geneva code until the red supergiant phase. 



-1 -0.5 0 0.5 1 

log (N/O) 



log (T.„) 


Fig. 10. Comparison between ON and non-ON fast rotating stars. Left: log(N/C) vs. log(N/0) diagram. Right: logg-log(r,.ff) diag ram. Large sym¬ 
bols are ON stars form the present study; open/small symbols are fa st rotating star s (V sin ; > 250 km s“‘) from iMartins et al.1 (20l5h. The startf Oph 
is also included: effective temperature and surface gravity are from iMartins et al] ( 120151) while surface abundances are from IVillamariz & Herrerd 
ll200.5h . 


spectroscopic binary (SBli lBolton & Rogerslll97^ . HD 14633 
an SBl with a low mass component and a high eccentricity 
(iBovaiian et al.ll20 05h. HD 10 2415 is reported as “likely vari- 
able by |Walborn et al.l (1201 ih and HD 201345 has an unclear 
status. 

Figure [TT] provides more information on the line-profile and 
radial-velocity variations of ten of the ON stars studied in the 
present paper. The data used to build Fig.[TT]are described in Ta¬ 
ble |AT| We chose Hei 45876 for study, as nearby telluric lines 
can be checked for wavelength-calibration issues. We confirm 
that HD 12323 and HD 14633 are SBl, with clear RV varia¬ 
tions. HD 48279 is very stable in our data, nor do HD 191423 or 
HD 201345 show signs of variability. All other stars examined 
(HD 13268, HD 91651, HD 102415, HD 117490, HD 123008) 


show some degree of RV variability, the origins of which are 
unclear (RV modulations, pulsations, wind variability). 

The exact position in respect of binarity among ON star is 
thus not completely clear. Variations in radial velocities are ob¬ 
served in a large fraction of the sample, but these are unambigu¬ 
ously attributable to binarity only in the few cases for which 
an orbital solution can be performed (HD 12323, HD 14633) 
or when an SB2 spectrum is observed (HD 89137). For the 
majority of the sample stars, radial velocity-variations do not 
have an unambiguous origin (they may be due to binarity or 
wind/photospheric variability), while a few stars’ spectra are 
consistent with constant radial velocities. 

This raises the question of the role of binarity in the strong 
chemical mixing that we have shown to be a clear-cut charac¬ 
teristic of ON stars. More specifically, one may wonder if mass 
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transfer is responsible for this strong enrichment. In the classi¬ 
cal mass-transfer scenario, material processed through nucleo¬ 
synthesis in the core of a companion is dumped onto what we 
observe today as an ON star. Stars with an SB 1 spectrum may be 
consistent with this scenario; the secondary may be a low-mass 
star, possibly a compact object resulting from the evolution and 
supernova explosion of an initially more massive star than the 
present-day primary. This initially more massive star could have 
contaminated the companion that we observe today as chemi¬ 
cally enriched before exploding as supernova. The system would 
have had to stay bound during the supernova phase to account 
for the SB 1 RV curve. A mass-transfer scenario may also be ap¬ 
plicable to RV-constant ON stars. These stars could have been 
initially part of binary systems in which mass transfer occurred 
before disruption in the supernova phase. 

iBolton & RogersI (Il978h argued that no sign of current 
episodes of mass transfer was observed among ON stars, espe¬ 
cially dwarfs. Plaskett’s star is a famous binary sys tem which is 
very p robably in a post Roche-lobe oyerflow phase (iLinder et alJ 
l2008h and which shows strong line emission and yariability. This 
is not obseryed currently among the ON stars we studied in the 
present paper, indicating that actiye mass transfer is probably 
not important among them. HD 48279, which has yery similar 
parameters to the SBl system HD 14633, does not show RV 
yariations. Both stars are located yery close to each other in the 
logg-Teff diagram (Fig. 0, hence their surface chemical enrich¬ 
ment cannot be clearly attributed to mass transfer, giyen that one 
is in a binary system and the other apparently is not, although we 
cannot exclude the possibility that HD 48279 is the result of the 
merger of two stars, which could explain, at least qualitatiyely, 
the strong mixing of CNO material. 

Howeyer, mass transfer is yery probably taki ng place in the 
binary BDh- 36 4063 (not included in our samnlel. [Williams et al.l 
(l2009l) showed that this system contains a ON9.7Ib star with 
periodic radial-yelocity yariations. The secondary component is 
not obseryed, although for its estimated mass it should be seen in 
high-resolution spectra. [Williams et al.l argue that the ON star is 
currently transferring mass to the companion, surrounding it by 
a thick disk, preyenting the secondary from being detected spec¬ 
troscopically. In this system, therefore, the ON star is currently 
the mass donor. 

A parallel situation is encountered in the massiye SB2 sys¬ 
tem LZ Cep. Unlike BDh- 36 4063, the companion in LZ Cep 
is clearly detected, excluding the presence of a th ick circum- 
stellar disk arising from ongoing mass transfer. iMahy et ^ 
(1201 ill obtained multi-epoch spectroscopy of the system and 
disentangled the spectra of both components, to which they as¬ 
signed spectral types of 09 III -i- ON9.7 V (primary, secondary). 
Their quantitatiye analysis of the indiyidual components’ spec¬ 
tra showed the secondary to be significantly nitrogen rich and 
carbon/oxygen poor: log(N/C) (respectiyely log(N/0)) reaches 
1.6 (1.3). The secondary component therefore shows the typi¬ 
cal ab undance patterns of ON stars, as established in the present 
study. iMahy et ^ (1201 ll) interpreted the surface abundances of 
the secondary as the result of a preyious episode of mass transfer 
in which the secondary was initially the most massiye star. The 
ON star in LZ Cep eyolyed faster and transferred mass to what 
was original the secondary (now seen as the primary). In this 
scenario, the abundance pattern of the ON secondary star is due 
to the remoyal of external layers during the mass-transfer phase, 
thereby exposing internal layers of the star. These internal layers 
are more chemically mixed and thus show CNO processed mate¬ 
rial; that is, for LZ Cep (and possibly BD-(-36 4063) the ON star 


would haye been created because of mass removal and not mass 
accretion. 

Thus, the role of binarity (and mass transfer) in the appear¬ 
ance of the ON phenomenon is not obyious. It is eyen less 
clear if we consider once again the fast rotator ^ Oph, which is 
a runaway (lBlaauwlll961l: iTetzlaff et al.ll20liT) . Runaways may 
be produced by dynamical interactions in young, dense clus¬ 
ters or by supernoya ejection from a binary system. Dynamical 
inter action in clusters usually inyolyes at least one binary sys - 

tem (iHut & BahcalJl983UHoffei<1983HHoogerwerf et al.l200lll . 

Consequently, there is a high probability that ^ Oph was part of 
a binary system at an ear lier stage i n its eyolution; in particular , 
iHoogerwerf et al.l (120001) (see also lyan Rensbergen et al.l [19961) 
argue that ^ Oph resulted from a supernoya kick. Whateyer the 
eyolutionary and dynamical history of this star, it did not reach 
the ley el of enrichment seen in ON stars. Consequently, the com¬ 
bination of binarity and fast rotation (possibly acquired through 
tidal spin-up) does not appear to lead ineyitably to the ON phe¬ 
nomenon. 

6. Conclusions and final remarks 

We haye performed a spectroscopic analysis of a sample of ON 
stars using atmosphere models computed with the code CMF- 
GEN. We haye determined the fundamental parameters and the 
He, C, N, and O surface abundances, along with projected rota¬ 
tional yelocities. The results can be summarized as follows : 

• These late-type ON stars haye initial masses in the range 
20 to 25 Mq, with the two supergiants of the sample being 
slightly more massiye. 

• The projected rotational yelocities of ON stars are on ayerage 
higher than those of comparison samples of norm al O stars 
(iMartins et al.ll2015t[^m(5n-DIaz & Flerre^l2014t) . 

• ON stars are chemically enriched; almost all haye a helium to 
hydrogen number ratio larger than 0.1. All stars are nitrogen 
rich, carbon poor, and, to a lesser extent, oxygen poor. The 
surface CNO abundances are consistent with nucleosynthesis 
predictions. 

• ON stars are more chemically mixed than morphologically 
normal stars of similar spectra l types and luminosity classes 
in the sample of iMartins et al.l (l2015l) . 

• Some ON stars are members of binary systems; some show 
radial-yelocity yariations of unclear origin; some are con¬ 
stant in radial yelocity. 

• Eyolutionary models including rotation are not able to re¬ 
produce the high degree of chemical mixing obseryed on the 
main sequence among ON stars. 

Erom these results, we conclude that ON stars are O stars 
showing CNO-processed material at their surface; the N/C and 
N/O ratios we establish are the highest obseryed so far in O stars. 
The processed material is the product of nucleosynthesis, but the 
mechanisms by which it is brought to the surface remains un¬ 
clear (rotational mixing and binary mass transfer being the prime 
suspects). ON stars rotate on ayerage faster than normal stars, but 
there exist fast rotators with similar masses and ages that are not 
ON stars. There are also binary and (presumably) single stars 
among the ON category, as well as non-ON fast rotators that are 
(probably) former members of binary systems. 
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Appendix A: Additional observational information 

Table lA.ll provides information on the spectra used to investi¬ 
gate spectral variability and binarity in the ON stars sample. The 
spectra are shown in Fig.ITTI 
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Fig. B.l. Best fit (red) of the observed spectrum (black) of HD 12323. 
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Fig. B.2. Best fit (red) of the observed spectrum (black) of HD 13268. 


Appendix B: Best fits to the observed spectra 
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Table A.l. Spectroscopic data used to study variability. 


Star 

date of observation 

Instrument 

Status 

HD 12323 

08 sep 2011 

N-NOT/FIES 

SBl 


12 sep 2011 

N-NOT/FIES 



29 oct 2012 

M-MERCATOR/HERMES 



29 oct 2012 

M-MERCATOR/HERMES 



25 dec 2012 

N-NOT/EIES 


HD 13268 

12jan 2011 

N-NOT/EIES 

variable 


26 oct 2012 

M-MERCATOR/HERMES 



28 jan 2013 

N-NOT/EIES 



29jan 2013 

N-NOT/EIES 


HD 14633 

13 jan 2009 

N-NOT/EIES 

SBl 


10 sep 2011 

N-NOT/EIES 



26 oct 2012 

M-MERCATOR/HERMES 



23 dec 2012 

N-NOT/EIES 


HD 48279 

14 jan 2011 

N-NOT/EIES 

constant 


24 dec 2012 

N-NOT/EIES 



25 dec 2012 

N-NOT/EIES 


HD 91651 

04 apr 2009 

ES02.2/EER0S 

variable 


11 feb2011 

ES02.2/EER0S 



20 mar 2011 

ES02.2/EER0S 



14 may 2011 

ES02.2/EER0S 



15 may 2011 

ES02.2/EER0S 



16 may 2011 

ES02.2/EER0S 


HD 102415 

12 may 2008 

ES02.2/EER0S 

variable 


16 may 2011 

ES02.2/EER0S 



17 may 2011 

ES02.2/EER0S 


HD 117490 

12 may 2008 

ES02.2/EER0S 

variable 


13 may 2011 

ES02.2/EER0S 



20 mar 2011 

ES02.2/EER0S 



21 mar 2011 

ES02.2/EER0S 



14 may 2011 

ES02.2/EER0S 



15 may 2011 

ES02.2/EER0S 


HD 123008 

13 may 2008 

ES02.2/EER0S 

variable 


16 may 2011 

ES02.2/EER0S 


HD 191423 

29aug 2011 

N-NOT/EIES 

constant 


11 sep 2011 

N-NOT/EIES 


HD 201345 

09 sep 2010 

N-NOT/EIES 

constant 


17jun 2011 

M-MERCATOR/HERMES 



10 sep 2011 

N-NOT/EIES 



25 dec 2012 

N-NOT/EIES 
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HD12323 HD13268 
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Fig. 11. Variability of ON stars. For each object, a selection of spectra centered on He 1 25876 is shown. The spectra are described in Table lA.fl 
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Fig. B.4. Best fit (red) of the observed spectrum (black) of HD 48279. 


Fig. B.5. Best fit (red) of the observed spectrum (black) of HD 91651. 
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Fig. B.6. Best fit (red) of the observed spectrum (black) of HD 102415. 
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Fig. B.7. Best fit (red) of the observed spectrum (black) of HD 117490. 


Fig. B.9. Best fit (red) of the observed spectrum (black) of HD 150574. 



Fig. B.8. Best fit (red) of the observed spectrum (black) of HD 123008. 


Fig. B.IO. Best fit (red) of the observed spectrum (black) of HD 191423. 
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Fig. B.ll. Best fit (red) of the observed spectrum (black) of HD 191781. 
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Fig. B.12. Best fit (red) of the observed spectrum (black) of HD 201345. 
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